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Abstract 

We report a novel method of determining the average Neel relaxation time and its temperature 
dependence by calculating derivatives of the measured time dependence of temperature for a frozen 
ferrofluid exposed to an alternating magnetic field. The ferrofluid, composed of dextran-coated 
FeaCR nanoparticles (diameter 13.7 nrn ±4.7 nm), was synthesized via wet chemical precipita¬ 
tion and characterized by x-ray diffraction and transmission electron microscopy. An alternating 
magnetic field of constant amplitude (Hq = 20 kA/m) driven at frequencies of 171 kHz, 232 kHz 
and 343 kHz was used to determine the temperature dependent magnetic energy absorption rate 
in the temperature range from 160 K to 210 K. We found that the specific absorption rate of the 
ferrofluid decreased monotonically with temperature over this range at the given frequencies. From 
these measured data, we determined the temperature dependence of the Neel relaxation time and 
estimate a room-temperature magnetocrystalline anisotropy constant of 40 kJ/m 3 , in agreement 
with previously published results. 
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INTRODUCTION 


Colloidal suspensions of superparamagnetic iron-oxide nanoparticles, specifically those of 
magnetite (FesCH) and maghemite ( 7 -Fe 2 C> 3 ), have been extensively investigated for their 
potential applications such as cell separation, use as contrast agents for magnetic resonance 
imaging, targeted drug delivery and magnetic fluid hyperthermia (MFH) 1, ^|. Among 
these applications, MFH has piqued the interest of researchers from various disciplines in¬ 
cluding biophysics, biomedicine and oncology due to its potential applications in various 
technologies for the treatment of cancer without the side effects inherent to radiation and 
chemotherapy-based methods 3]. MFH involves the excitation of magnetic nanoparticles 
suspended in a fluid medium (a ferrofluid) using an oscillating magnetic field of the form 
H(t) = ihocos(27r/t), where H 0 is the field amplitude (typically between 5-30 kA/m) and 
/ is the frequency (typically in the range from 150 kHz - 350 kHz). For a single-domain 
nanoparticle, hysteresis losses are absent and energy absorption from the field can occur 
via two excitation mechanisms: Neel and Brownian relaxation. This absorbed magnetic en¬ 
ergy is eventually transformed to thermal energy and the temperature of the ferrofluid rises. 
Neel relaxation [4] involves the alignment of the nanoparticle moment with the external field 
within a fixed, non-rotating nanoparticlc. As the excitation occurs against the anisotropy 
energy barrier of the particle the relaxation time depends strongly on the nanoparticle’s 
magnetic volume, V m , magnetocrystalline anisotropy constant, K , temperature, T and char¬ 
acteristic relaxation time, tq. While there are a few models to describe this relaxation, the 
most commonly used has the form [4] 
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Brownian relaxation [5[ involves the alignment of the particle’s moment with the external 
field via the physical rotation of a fixed-moment nanoparticle within the carrier fluid. This 
alignment is affected by the hydrodynamic properties of the carrier fluid and nanoparticles, 
it is described by [6j| 
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where r] is the viscosity of the carrier fluid and V B is the hydrodynamic volume of the 
composite particle which includes any surfactant layer used to give colloidal stability to 
the ferrofluid. When both mechanisms are active these processes occur in parallel and the 
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effective relaxation time which describes the energy transfer rate is given by 

t n t b 
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From this relation, it is evident that the shorter relaxation time dominates the overall 
dissipative characteristics and thus determines the heating rate of the sample. 

The Neel mechanism plays a dominant role in the relaxational characteristics of the 
nanoparticles particularly when they are embedded in a cancerous tissue undergoing MFH 
treatment, since in the tissue’s environment, due to immobilization of the particles, the 
Brownian mechanism is highly damped. Several studies have shown that if the magnetic 
nanoparticles are internalized by the cancer cells, they are either locked in the cell plasma or 
adhere to the cell walls [7]. Consequently, the role of Brownian relaxation is insignificant for 
nanoparticles used in the hyperthermia treatment of cancer when the ferrofluid is applied 
directly to the tumor tissues via injection. In such situations, Neel relaxation produces 
the local heating of the cancerous mass; however, one must be careful about comparing 
experiments done in the laboratory with therapies that will be performed under real-life 
conditions. 

The experimental studies of the heating rates of ferrofluids in the presence of alternating 
magnetic fields can be broadly placed into one of two categories: those focusing on varying 
only the intrinsic physical and magnetic parameters of the nanoparticles 8, 9] and those 
that investigate the influence of varying only the parameters of the applied field (i.e. H 0 
and /) Q. Among these two broad classifications, studies of how the intrinsic parameters 
of nanoparticles influence magnetic heating still remain an active area of research. However, 
as the effects of the intrinsic parameters are reflected through the relaxation processes (Neel 
and Brownian) which act in parallel, it is difficult to understand the contribution of each 
process on the heating characteristics of the ferrofluid. 

Knowing the significance of the Neel mechanism in hyperthermia, many studies have been 
performed to estimate this parameter using ac and dc magnetic susceptibility measurements 
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12]; however, these methods often require large volumes of sample, accurate determi¬ 


nations of sample volume and mass and costly equipment. In this paper, we present a new 
method to determine this parameter using a simple induction heating system. This method 
requires only measurements of the sample heating rate under the influence of an oscillating 
magnetic field of constant amplitude at two different frequencies. Additional measurements 
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of sample mass, volume and magnetic anisotropy are not needed. The goal of this investiga¬ 
tion is to determine the temperature dependence of the Neel relaxation time (tn) when the 
Brownian mechanism is completely quenched. This suppression of the Brownian relaxation 
is achieved by performing the experiments in the temperature range from 160 K to 210 K 
where the carrier fluid (DI water) is completely frozen, thus locking the nanoparticles in 
place. 

Experimentally, the heating rate in MFH is expressed in terms of the specific absorption 
rate ( SAR ) which is defined to be the power absorbed per unit mass of the nanoparticles 
in the ferrofluid. This can be calculated using the thermodynamic relation for rate of heat 
energy absorbed by the sample as [ip| 
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where M samp i e is the overall mass of the sample, mFe 3 o 4 is the mass of the Fe 3 0 4 nanoparti¬ 
cles, C(T ) is the temperature dependent specific heat of the carrier fluid and ATheating /At 
is the time rate of change of the sample’s temperature as it absorbs energy from the applied 
magnetic field. When the sample begins to warm under the influence of the applied mag¬ 
netic field, it is important to note that heat exchange with the environment is also occurring 
due to the temperature differential between the sample and its surroundings. For accurate 
determination of SAR, it is important to estimate this heat exchange between the sample 
and the environment due to convective, conductive and radiative processes. While this pro¬ 
cess is quite hard to do from a calculational standpoint, we can experimentally estimate 
the collective effects of heat exchange via these processes by measuring the heating of the 
sample placed into the experimental apparatus while keeping the external field turned off. 
From this data, we determine the specific power gain (SPG) as heat is transferred from the 
environment to the sample as [t3 ] 
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/At is the warming rate of the sample due only to heat exchange between 


the sample and the environment. Once this SPG is determined, we can correct the measured 
SAR for the sample as a function of temperature jpll ] 
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Due to the suppression of the Brownian mechanism, the values of SAR obtained over this 
temperature interval are dependent only upon Neel relaxation. This relaxation time can 
be determined using a theoretical model based on a direct relationship between the power 
density dissipated by the nanoparticles and the out-of-phase, dissipative component of the 
ferrofluid’s susceptibility [15 ] 

2nfT N 


P = nfJoHoXof- 


(7) 


1 + (27r/Tjv) 2 ’ 

where Xo is the equilibrium susceptibility and fi 0 is the permeability constant. In this model, 
the energy dissipation by the nanoparticles, which is related to the out-of-phase component 
of the magnetic susceptibility, is expressed in terms of the equilibrium susceptibility and the 
Neel relaxation time tm- Measuring the SAR at two different frequencies f\ and f 2 while 
keeping the field amplitude (H 0 ) fixed, we find Tn from Equation [7] as 
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It is interesting to note that tjv is determined only from the values of / 1; f 2 and the sam¬ 
ple heating rates (corrected for heat exchange with the environment) and do not require 
measurement of the sample mass, nanoparticle mass, or specific heat values. 


EXPERIMENTAL DETAILS 

Iron oxide nanoparticles were synthesized using a standard co-precipitation technique in 
which an aqueous solution of FeCl 3 - 6 H 2 0 and FeCl 2 -4H 2 0 were mixed in a 2:1 molar ratio 
and Fe 3 C >4 nanoparticles were precipitated by the drop-wise addition of 1M NH 4 OH. During 
precipitation, N 2 gas was bubbled through the solution to protect against oxidation of the 
Fe 2+ ions into Fe 3+ ions. The precipitate was separated from the solution by a strong magnet, 
washed with DI water and re-suspended in a metastable 0.5 M NaOH solution. In order to 
suspend the precipitated nanoparticles in a carrier solution (DI water) they were coated in 
dextran by the drop-by-drop addition of the metastable solution of Fe 3 C >4 to a solution of 
15-20 kDa dextran (MP Biomedicals) in 0.5M NaOH while simultaneously probe sonicating. 
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The product was rinsed and resulted in a water-based suspension of Fe 3 04 nanoparticles 
with a concentration of 20 mg of Fe 3 C >4 per mL of solution. A portion of the sample was 
lypholized and characterized via x-ray diffraction (Rigaku MiniFlex 600) and transmission 
electron microscopy (JEOL HR TEM 2010 operating at 200 keV). 

For calorimetric measurements in the 160 K to 210 K range, the ferrofluid sample was 
first cooled to 77 K via the immersion of seald vials in liquid nitrogen. The sample was 
then allowed to warm under ambient conditions to 160 K at which time a 20 kA/rn (250 
Oe) ac magnetic field was applied using an Ambrcll EasyHeat 2.4 kW induction heating 
system with a water-cooled, 8-turn, 2-cm-diameter coil. Temperature versus time data 
were collected in the 160 K to 210 K region using an Optocon FOTEMP1-H fiber optic 
temperature monitoring system equipped with a TS5 optical temperature sensor with 0.1 K 
accuracy. The sample vial was thermally insulated using cotton padding and foam rubber 
to minimize heat exchange with the environment; however, as thermal interaction with the 
environment cannot be completely suppressed, experiments were performed to determine 
the extent of the heat exchange with the surrounding environment and this estimation was 
accounted for in all measurements as discussed above. 

RESULTS AND DISCUSSION 

The left panel of Figure |T] shows the powder x-ray diffraction (XRD) spectrum taken 
from the lypholized dextran coated Fe 3 04 nanoparticles. The open symbols represent the 
observed counts recorded for different d-spacing values between 1.25 and 3.25 A. The solid 
line is a full-profile Le Bail fit to the data, the vertical bars indicated the d-spacing positions 
of the Bragg reflections and the lower trace is the difference curve between the observed and 
calculated counts. The fit confirms that the sample consists of a single nanocrystallinc phase 
of Fe 3 C >4 with cubic FmZm symmetry and lattice constant a = 8.36 A. Using the full-width 
at half maximum of the (311) reflection in Scherrer’s equation, an average nanoparticle 
diameter of 14 nm was estimated for this sample. The right panel of Figure [Tj shows the 
magnetic field dependence of the ferrofluid’s magnetization measured at 150 K which exhibits 
no hysteresis thus confirming the superparamagnetic nature of the nanoparticles within the 
experimental temperature range between 160 K and 210 K. 

Transmission electron microscopy (TEM) was performed to determine the mean diameter 


6 


and size distribution of the nanoparticles in addition to confirming the crystallographic 
information determined via XRD. Figure [2](a) shows a histogram of particle sizes determined 
from the bright-held TEM micrograph shown in Figure [2Xb) and similar images at the same 
magnification. The histogram shows a log-normal particle size distribution having mean 
diameter ( D) = 13.4 nm and standard deviation from the mean of an = 4.7 nm. This number 
is in agreement with the 14 nm average diameter determined using Scherrer’s equation in 
conjunction with the collected XRD spectrum. Figure [2](c) shows a high-resolution TEM 
micrograph of a portion of a single nanoparticle. The visible lattice planes in the image 
show a spacing of 2.9 A marking them as the (220) set of planes in FesCR. Lastly, Figure 
[2(d) shows a selected area electron diffraction (SAED) pattern which was used to to confirm 
the cubic symmetry determined using XRD. 

The left panel of Figure [3] shows the temperature versus time data collected while heating 
the ferrofluid using an alternating magnetic held of amplitude H 0 = 20 kA/m at 171 kHz 
(open circles), 232 kHz (open triangles) and 343 kHz (open squares). All three curves 
were fit to a polynomial (solid lines) and differentiated in order to determine the slope, 
AT/At. To account for any heat exchange with the environment, a control experiment was 
performed in which the ferrofluid was placed inside the held coil and allowed to warm in the 
absence of the magnetic held. This data is plotted in the right panel of Figure [3] and will 
be referred to as the ambient curve. From this data, we were able to calculate the specihc 
absorption rate ( SARh ea ung ) for each frequency and the specihc power gain ( SPG warm i ng ) 
for the ambient experiment all as functions of temperature using the temperature-dependent 
specihc heat of ice 16] and Equations [4] and [5] respectively. From these calculations, we 
were able to determine the specihc absorption rate (SAR) due only to energy absorbed from 
the alternating magnetic held using Equation [6l The results of these calculations are shown 
in Figure [4] and yield average values of SAR over the temperature interval of 57 W/g, 88 
W/g and 130 W/g at 171 kHz, 232 kHz and 343 kHz, respectively. 

Additionally, using the temperature dependent slopes of the data presented in Figure [3] 
found through differentiation of the polynomial hts, we are able to use Equations El and [9] 
to calculate the temperature dependence of the Neel relaxation time over the temperature 
regime in question. The results of this calculation are shown in Figure [5] and show a value 
near 5 x 10” 7 s for all possible permutations of f\ and fi- Using Equation [Tj with this 
determined value of tjv, the experimentally determined particle diameter of 13.7 nm, and 
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assuming a characteristic relaxation time of To = 1CU 9 s we find the room temperature 
magnetocrystalline anisotropy constant, K , for this sample to be 40 kJ/m 3 , which agrees 
with the range of values previously reported in the literature HQ. 

It is important to emphasize that for quantitative determination of tjv using Equation |H|, 
one needs the time derivatives of the temperature at two different frequencies and for the 
ambient (no applied field) case. Using these values in Equation [9l we determine the temper¬ 
ature dependence of the Neel relaxation time using Equation El In our case, by taking data 
at three different frequencies, we were able to provide three different estimations of the tem¬ 
perature dependence of tjv and found, as expected, that the value was roughly the same for 
each pair of frequencies. We see in Figure 0 a convergence of the values of r N as we approach 
210 K. In addition, it is important to note that the magnitude of the Neel relaxation time 
is on the order of 10 -7 s as expected for a system of non-interacting magnetic nanoparticles 
having characteristic time constant To~lCU 9 — 1CU 13 s [l, 19]. The non-interacting nature 
of this ensemble was conhrmed through calculation of the relative variation of the blocking 
temperature per frequency decade (0 = AT/Tlog 10 /) in which we found a value of </> = 0.18, 
an order of magnitude larger than that found in interacting samples exhibiting spin-glass-likc 


transitions 


20H23j]. The details of this investigation will be published in a subsequent paper 


24], At the macroscopic level, tn, as described by Equation [T] references the relaxation of 


an individual nanoparticle of volume V m having magnetocrystalline anisotropy constant K 
and characteristic relaxation time To = \JnKV m J4k bT ; however, in a polydisperse system 
these parameters may vary from particle to particle and consequently their Neel relaxation 
times do vary. In SAR measurements, one measures the temperature rise of the ferrofluid 
due to the average thermal energy of its constituents (i.e. nanoparticles, surfactant and 
carrier fluid) in contrast to the relaxation time described by Equation [T] which refers to a 
single particle. Thus, our measured represents an average Neel relaxation time of the 
ensemble of particles-this is what is responsible for the measured values of SAR. However, 
for a monodisperse system, Equations [j] and El may be applied as what is true for any one 
particle is true for any other in the system. The use of Equation [U however, requires an 
accurate determination of the particle volume, V m , and the magnetocrystalline anisotropy 
constant, K. On the other hand, the method outlined in this paper requires only the value 
of the applied field frequencies and the derivatives of the temperature versus time curves 
which can be accurately determined using simple to understand methodology. 












CONCLUSION 


We studied the temperature dependence of the specific absorption rate ( SAR ) of a frozen 
ferrofluid at different frequencies at a constant field amplitude. From the determined SAR 
values we were able to calculate the temperature dependence of the Neel relaxation time 
using a novel and simple approach. This method is useful for cases in which there are 
uncertainties related to the magnetocrystalline anisotropy constant, K, and nanoparticle 
size. In our experimental temperature window between 160 K and 210 K, our measured 
SAR decreases with temperature; however, this trend may not continue as the experimental 
window is widened to include temperatures above 210 K where other effects may begin 
to play a role in the relaxation characteristics of the nanoparticles. We believe that this 
investigation offers a new and simpler method of determining an average Neel relaxation 
time in a colloidal suspension of magnetic nanoparticles. 
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FIG. 1: (Left Panel) A powder x-ray diffraction spectrum of a lypholized portion of the sample. 
The open symbols correspond to the oberserved data and the solid line is a full-profile Le Bail 
fit to the data indicated the Fm2>m cubic symmetry and 8.36 A lattice constant characteristic 
of Fe 3 C> 4 . The Miller indices of the more intense reflections are indicated and the bottom trace 
is the difference curve between the observed data and the fit. (Right Panel) The magnetic field 
dependence of the sample’s magnetization at 150 K showing that the nanoparticles are in the 
superparamagnetic state in the experimental temperature regime between 160 K and 210 K. 
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FIG. 2: (a) A histogram of the particle size distribution determined using the micrograph shown 
in (b) and other similar micrographs taken at the same magnification, (c) A high-resolution TEM 
micrograph of a portion of a single nanoparticle with an outline of the particle provided for clarity. 
The spacing of the observed lattice planes was found to be 2.9 A indicative of the (220) set of 
planes, (d) A selected area electron diffraction (SAED) pattern which was used to confirm the 
cubic nature of the sample. 
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FIG. 3: (Left Panel) Temperature versus time measurements made with the ferrofluid in an alter¬ 
nating magnetic field of amplitude Hq = 20 kA/rn at frequencies of 171 kHz (A), 232 kHz (B) and 
343 kHz (C). (Right Panel) Temperature versus time measurements made with the ferrofluid in 
the field coil but in the absence of an applied field. 
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FIG. 4: The temperature dependence of the specific absorption rate (SAR) of the ferrofluid mea¬ 
sured at frequencies of 171 kHz (open circles), 232 kHz (open triangles) and 343 kHz (open squares). 
The temperature dependence of the specific heat of ice used to calculate the specific absorption 
rate is plotted as a solid line 
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FIG. 5: The temperature dependence of the Neel relaxation time calculated using equations 8 and 
9 for frequencies /i = 171 kHz and /2 = 232 kHz (open circles), /i = 171 kHz and /2 = 343 kHz 
(open triangles) and /i = 232 kHz and = 343 kHz (open stars). 
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